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1. Introduction 

The main task in application of the quantum mechanics is to solve the Schrodinger 
equations with different potentials. Unfortunately, realistic physical problems can 
practically never be solved exactly. Then one has to resort to some approximations. 
Most widely used among them is the perturbation theory. However, the explicit 
calculation with the Rayleigh - Schrodinger perturbation theory, described in most 
quantum mechanics textbooks, runs into the difficulty of the summation over all 
intermediate unperturbed eigenstates. To avoid this difficulty, alternative perturbation 
procedures have been proposed. They are Sternheimer's method [1,2], the Dalgarno- 
Lewis technique [3-6], the method developed by Zel'dovich [7-9] and the logarithmic 
perturbation theory [10-17]. 

It must be stressed that in both the Steinheimer method and the Dalgarno- 
Lewis technique one has to solve inhomogeneus differential equations. In solving the 
Schrodinger equation with the Zel'dovich method, for avoiding the summation over 
intermediate states the Lagrange condition of the theory of differential equations is 
used. 

The logarithmic perturbation theory is more straightforward in the sense that it 
does not require solving any equation or involving any additional condition. Within 
the framework of this approach, the conventional way to solve a quantum-mechanical 
bound-state problem consists in changing from the wave function to its logarithmic 
derivative and converting the time-independent Schrodinger equation into the nonlinear 
Riccati equation. 

In the case of ground states, the consequent expansion in a small parameter 
leads to handy recursion relations that permit us to derive easily the corrections to 
the energy as well as to the wave function for each order. Notice that the obtained 
series are typically divergent and the evaluation of perturbation terms of large orders is 
needed for applying the modern procedures of summation of divergent series. However, 
when radially excited states are considered, the standard technique of the logarithmic 
perturbation theory becomes extremely cumbersome and, practically, inapplicable for 
describing higher orders of expansions. 

For this reason authors of the paper [9] conclude that in dealing with excited states 
the method developed by Zel'dovich has a remarkable advantage over the logarithmic 
perturbation theory. But in fact there is another approach to obtaining expansions of 
the logarithmic perturbation theory that describes excited states and the ground state 
exactly in the same manner by means of simple recursion formulae. 

Indeed, the above mentioned disadvantage of the standard approach to the 
logarithmic perturbation theory is caused by factoring out zeros of the unperturbed 
wave functions with taking into account corrections to the positions of these nodes. 
On the other hand, the number of zeros of the wave function most conveniently and 
naturally is introduced in the consideration by means of quantization conditions for the 
logarithmic derivative of the wave function just as within the framework of the Wentzel- 
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Kramers-Brillouin (WKB) method [18-20]. However, since the WKB-approximation is 
more suitable for obtaining energy eigenvalues in the limiting case of large quantum 
numbers but the perturbation theory, on the contrary, deals with low-lying levels, the 
WKB quantization conditions need change. 

Recently, a new procedure based on specific quantization conditions has been 
proposed to get series of the logarithmic perturbation theory via the ^-expansion 
technique within the framework of the one-dimensional Schrodinger equation [21]. 
Avoiding disadvantage of the standard approach, this straightforward semiclassical 
procedure results in new handy recursion formulae with the same simple form both 
for the ground state and excited states. Moreover, these formulae can be easily applied 
to any renormalization scheme of improving the convergence of expansions [22]. 

However, in most of the practical applications of quantum mechanics, one deals 
with the more complicated case involved the three-dimensional Schrodinger equation 
with the anharmonic oscillator potential or the screened Coulomb potential. The object 
of this paper is to extend the above mentioned formalism to the bound-state problem 
for the spherical anharmonic oscillator that has numerous applications in the theory 
of molecules and solid state physics. The another widely used bound-state problem 
for the three-dimensional Schrodinger equation with the potential, having the Coulomb 
singularity, will be published elsewhere. 

The layout of the paper is as follows: in the next section we summarize the main 
ideas of our approach concerning the construction of new quantization conditions which 
are more suitable for the semiclassical treatment of the logarithmic perturbation theory. 
In section 3 the quantization conditions obtained are used for deriving new simple 
recursion relations for the calculation of perturbation expansions, which have the same 
simple form both for the ground state and excited states. Section 4 demonstrates that 
in the case of the harmonic oscillator potential the described approach restores the exact 
results for the energy eigenvalues and eigenfunctions. Here we illustrate our method by 
applying it to the example of the quartic anharmonic oscillator, as well. The paper ends 
with a brief summary. 

2. The classical limit and the quantization rule 

We study the bound-state problem for a non-relativistic particle moving in a central 
potential of an anharmonic oscillator admitted bounded eigenfunctions and having in 
consequence a discrete energy spectrum. This potential has a single simple minimum 
at the origin and is given by a symmetric function V(r) which can be written as 

V(r) = -mu 2 r 2 + Y J V^r■ 2l+2 ■ (1) 
2 i>i 

Then, by separating the angular part, the reduced radial part of the Schrodinger 
equation takes the form 

-^W+( ? l^ + w)l/(r) = «/(r). (2) 



Logarithmic perturbation theory for the spherical anharmonic oscillator 



4 



Following usual practice, we apply the substitution, C{r) = hU'(r)/U(r), accepted 
in the logarithmic perturbation theory and go over from the Schrodinger equation (J2J) 
to the Riccati equation 

hC'(r) + C 2 {r) = - Z ^ 2 + ^ + 2mV{r) - 2mE. (3) 

We attempt to solve it explicitly in a semiclassical manner with series expansions 
in the Planck constant 

oo oo 

E = J2E k h k C(r) = J2C k (r)h k . (4) 

The ^-expansions under discussion simplify the problem of taking into account the nodes 
of wave functions for excited states, allowing the use of the quantization condition and 
the formalism of the theory of functions of a complex variable. 

In the complex plane, a number of zeros iV of a wave function inside the closed 
contour is defined by the principle of argument known from the analysis of complex 
variables. Being applied to the logarithmic derivative, C(r), it means that 




This quantization condition is exact and is widely used for deriving higher-order 
corrections to the WKB-approximation [23,24] and the 1/iV-expansions [25-27]. There 
is, however, one important point to note. Because the radial and orbital quantum 
numbers, n and I, correspondingly, are specific quantum notions, the quantization 
condition (JHJ must be supplemented with a rule of achieving a classical limit for these 
quantities. It is this rule that stipulates the kind of the semiclassical approximation. 

In particular, within the framework of the WKB-approach the passage to the 
classical limit is implemented using the rule 

H — > n — > oo I — > oo hn = const hi = const (6) 

whereas the 1/iV-expansion requires the condition [25-27] 

h~ — ► n = const I — > oo hn —>■ hi = const. (7) 

The semiclassical treatment of the logarithmic perturbation theory proved to 
involve the alternative possibility: 

h~ — ► n — const I = const hn — > hi — *■ 0. (8) 

Let us consider the latter rule from the physical point of view. Since hi — > as 
% — ► 0, the centrifugal term, h 2 l (I + 1) /r 2 , has the second order in h and disappears 
in the classical limit that corresponds to falling a particle into the center. This means 
that a particle drops into the bottom of the potential well as H — *■ and its classical 
energy becomes Eq = minV r (r) = 0. It appears from this that the series expansion in 
the Planck constant for the energy eigenvalues must now read as E = J2k^=i Ekh k . 
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Upon inserting the ^-expansions for E and C(r) into the Riccati equation J2J) and 
collecting coefficients of equal powers of h, we obtain the following hierarchy of equations 

Clir) = 2 mV(r) 
C' (r)+2C (r)C 1 (r) = -2 mE 1 

C[{r) + 2 C (r)C 2 (r) + C^r) = ^±22 _ 2 m £ 2 (9) 

k 

C'k-i(r) + E CiCOCfc-iCr) = -2 m£ fc fc > 2. 

i=0 

In the case of ground states, the recurrence system at hand coincides with one 
derived by means of the standard technique and can be solved straightforwardly. For 
excited states, however, it is necessary to take into account the nodes of the wave 
function, that we intend to do by making use of the quantization condition (jSJ). 

It should be stressed that our approach is quite distinguished from the WKB 
method not only in the rule of achieving a classical limit but also in the choice of 
a contour of integration in the complex plane. With a view to elucidating the last 
difference let us now sketch out the WKB treatment of the bound-state problem for the 
case of the spherical anharmonic oscillator. In the complex plane, because the potential 
is described by the symmetric function ((H), this problem has two pairs of turning points, 
i.e. zeros of the classical momentum. Therefore we have two cuts between these points: 
in the region r > as well as in the region r < 0. In spite of only one cut lies in the 
physical region r > 0, the contour of integration in the WKB quantization condition has 
to encircle both cuts for the correct result for the harmonic oscillator to be obtained [28]. 

In our approach, when a particle is dropping into the bottom of the potential well 
these four turning points are drawing nearer and, at last, are joining together at the 
origin. Hence, all nodes of the wave function are now removed from both positive and 
negative sides of the real axis into the origin and our contour of integration must enclose 
only this point and no other singularities. 

Further, let us count the multiplicity of a zero formed at r = 0. For the regular 
solution of the equation (J2J), the behaviour r l+1 as r —> brings the value I + 1. The 
number of nodes of eigenfunction in the region r > is equal to the radial quantum 
number n. But because the potential (0) is a symmetric function the same number 
of zeros must be in the region r < 0, too. Then the total number of zeros inside the 
contour becomes equal to iV = 2n + I + 1. 

Taking into account the first order in H of the right-hand side, the quantization 
condition (jHJ) is now rewritten as 

— <fd(r)dr = 2n + / + l — I Gdr) dr = jfe ^ 1. (10) 

27T1 J 27T1 J 

A subsequent application of the theorem of residues to the explicit form of functions 
Cfc(r) easily solves the problem of the description of the radially excited states. 
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3. Recursion formulae 

Let us consider the system and investigate the behaviour of the function Ck(r). From 
the first equation it is seen that Co(r) can be written in the form 

(q oo \ 1/2 oo 

i + — ^E^ r2i = r Y, c i r2t ( n ) 
muJ i=l / i=0 

where the minus sign is chosen from the boundary conditions and coefficients C? are 
defined by parameters of the potential through the relations 

C° = -mw C? = J- fa C° p Cl p - 2mv\ i > 1. (12) 

From (fTTj) we recognize that C o (0) = and, consequently, the function Ci(r) has a 
simple pole at the origin, while Ck(r) has a pole of order (2k — 1). Thus Ck{r) can be 
represented by the Laurent series 



CO 



C7*=- ' 



(3 - 2fc + 2^)^" 1 + EE C?Cj? + 2 E C?Ct P - Z(Z + l)6 2 , k 6 Q 



.(15) 



C fc (r) = r 1 - 2 ' 2 C-r 2i k > 1. (13) 

i=0 

With definition of residues, this expansion permits us to express the quantization 
condition (|10|) explicitly in terms of the coefficients C\ as 

Cti = (2n + Z + 1)5 M . (14) 

It is this quantization condition that makes possible the common consideration of 
the ground and excited states and permits us to derive the simple recursion formulae. 

The substitution of the series (JI2J) and (fT3j) into the system © in the case i ^ k — 1 
yields the recursion relation for obtaining the Laurent-series coefficients 

fc-l i 

^/-liz-ih-j , oV^^fO^yfc i/7 , i\r r 

0,i 

j=lp=0 p=l 

Hi = k — 1, by equating the explicit expression for G\_ x to the quantization 
condition (fT^j) we arrive at the recursion formulae defined the perturbation corrections 
to the energy eigenvalues 

2mE k = -Ct\ - E E qJCtt, • (16) 

3=0 p=0 

Thus, the problem of obtaining the energy eigenvalues and eigenfunctions for the 
bound-state problem for the anharmonic oscillator can be considered solved. The 
equations (JT5|l and ()16|) have the same simple form both for the ground and excited 
states and define a useful procedure of the successive calculation of higher orders of the 
logarithmic perturbation theory. 
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On applying the recursion relations obtained, the analytical expressions for first 
corrections to the energy eigenvalues of the spherical anharmonic oscillator ([T} are found 
to be equal to 
1 + 2N 

Ei = —2— w 

„ _ (3-2L + 6r/) Vl 

2 ~~ A 2 — 2 

4m z oo z 

E 3 = 1 + , 2J !" ((-21 + 9L- 1777) V + m (15-6L + 10 77) cu 2 v 2 ) 
£4 = 77; -a o f (333 + llL 2 -3L (67 + 86 77) + 3 77 (347+ 125 7?)) Vl 3 

-6m (60 + 3 (-13 + L) L + 175 77 — 42 L 77 + 55 t? 2 ) a; 2 1^2 (17) 
+m 2 (6 L 2 - 12 L (6 + 5 77) + 35 (3 + 2 77 (4 + 77))) a/ 1> 3 ) 

E 5 = 1 + 2 g iV 11 f f 30885 + 909 L 2 — 27 L (613 + 330 77) + 77 (49927 + 10689 77)) v x A 

128 777 s to 11 V v ' 

- 4 771 (11220 + 393 L 2 - 6 L (1011 + 475 77) + 77 (16342 + 3129 77)) uj 2 ?Ji 2 v 2 
+ 16 777 2 (33 L 2 — L (501 + 190 77) + 63 (15 + 77 (19 + 3 77))) uj 4 Vl v 3 
+ 2 m 2 (3495 + 138 L 2 + 4538 77 + 786 if - 30 L (63 + 26 77)) uj 4 v 2 2 

-4m 3 (30L 2 -20L (24 + 777) + 63 (15 + 277 (8 + 77))) u 6 v 4 

where N = 2 n + 1 + 1, 77 = N (N + 1), L = 1(1 + 1). 

As it is seen, the obtained expansion is the expansion in powers of the Taylor-series 
coefficients for the potential function. 

It is also evident that for the energy eigenvalues, when k = 1, we readily have the 
oscillator approximation [29] 

E 1 =(2n + l + I) oj. (18) 

Let us demonstrate that in the case of the isotropic harmonic oscillator our 
technique restores the exact solution for the wave functions, U(r), too. 

Putting, for simplicity, h = m = u = 1, from equations (|T2"|) . (fTB'j) and (fTKj) we find 

that 

C (r) = -r C k (r) = d k r l ~ 2k k>0 (19) 
where coefficients d k obey the relations 
do = -1 
di =N 

2 d 2 = N 2 - N - 1(1 + 1) (20) 



24 = 



fc-i 

(3 - 2k) 4_i + djd k _j k > 2 
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Due to the definition C(r) = U'(r)/U(r), the straightforward integration of the 
function Co(r) and the part of the function C\{r) immediately gives factors e _r2//2 
and r l+l in the eigenfunctions, U(r), providing their correct behaviour at infinity and 
the origin. The remaining part is a polynomial P n (r) = J2k=oPk r2k that satisfies the 
equation 

oo 

Pn(r 2 )/Pn(r 2 ) = £ - (/ + IK 1 . (21) 

k=l 

The polynomial coefficients, pk, are determined by the system 

n 

2 p m {n - m) + £ Pi <V m +i = ■ ( 22 ) 

The combination of these equations, multiplied by a suitable dj with a view to taking 
into account equation pUjl. arrives at the following relation between two consecutive 
coefficients: 

Pm-i = m(m + l + 1/2) ^ 23 ^ 
p m m — n — 1 

that is the known recursion formula for the associated Laguerre polynomial, L^ +1 / 2 (r 2 ), 
of degree n [30]. Thus, the described technique restores the exact result for the 
unnormalized wavefunctions of the harmonic oscillator [29], too. 

As an example of application we examine excited states as well as the ground states 
of the quartic anharmonic oscillator with a potential 

V(r) = mu 2 r 2 /2 + Ar 4 (24) 

where A is a positive constant. 

Equations (fTTjl for the energy eigenvalues now are rewritten as 

(3-2L + 6??) 

^2 = ~. — 5 A 

-(1 + 2AQ (21-9L + 17„) A2 

8 m 4 cij 5 

(333 + llL 2 -3L (67 + 86 77) + 3 1] (347+ 125 77)) 3 

= 77 — r — r ^ 

16 m° c<j 8 

-(1 + 2 iV) (30885 + 909 L 2 — 21 L (613 + 330 77) + 77 (49927+ 106897/)) ^ 4 
5 ~ 128 m s u n 

We recall that here N = 2 n + I + 1, 77 = JV (N + 1), L = 1(1 + 1). 

It is readily seen that the use of the ^-expansion technique does lead to the explicit 
perturbation expansion in powers of the small parameter A. 

In the case of ground states, obtained expansions for the energy eigenvalues coincide 
with those listed in [31]. In the case of excited states, our corrections coincide with 
corrections up to the second order which are just only calculated in [32]. 

Thus, we derive the recursion formulae for obtaining arbitrary order terms of the 
weak coupling power series for the bound-state problem of the anharmonic oscillator. 
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Unfortunately, as it is known, this series is asymptotic and diverges for any finite value of 
the parameter A that requires the use of some procedures of improving the convergence 
(for references see [33]). The most common among them are various versions of the 
renormalization technique, intended to reorganize a given series into another one with 
better convergence properties. In practice, these procedures involve quite a number 
of the perturbation terms and are stopped at a finite order approximation when the 
change of some eigenvalue with increasing order of approximation becomes less than the 
needed exactness. If the series begins diverge before achieving the needed exactness of 
the approximation, we must apply one of the methods of summation of divergent series. 

It should be noted, that the proposed formalism of the logarithmic perturbation 
theory is easily adapted to the treatment of any renormalization scheme in terms 
of handy recursion relations within the framework of the united approach [22]. 
Besides, the described technique can be extended to the case of the screened Coulomb 
potential through the modification of the quantization conditions, that will be published 
elsewhere. 

5. Summary 

In conclusion, a new useful technique for deriving results of the logarithmic perturbation 
theory has been developed. Based upon the /i-expansions and suitable quantization 
conditions, new handy recursion relations for solving the bound-state problem for a 
spherical anharmonic oscillator have been obtained. These relations can be applied to 
excited states exactly in the same manner as to ground states providing, in principle, the 
calculation of the perturbation corrections of large orders in the analytic or numerical 
form. Besides this remarkable advantage over the standard approach to the logarithmic 
perturbation theory, our method does not imply knowledge of the exact solution for zero 
approximation, which is obtained automatically. And at last, the recursion formulae at 
hand, having the same simple form both for the ground state and excited states, can be 
easily adapted to applying any renormalization scheme for improving the convergence 
of obtained series, as it is described in [22]. 

The another widely used case of the three-dimensional Schrodinger equation with 
the potential having the Coulomb singularity will be published elsewhere. 
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